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We analyse the full counting statistics (FCS) of the charge transport through the Anderson 
impurity model (AIM) and similar systems with a single conducting channel. The object of principal 
interest is the generating function for the cumulants of charge current distribution. We derive an 
exact analytic formula relating the FCS generating function to the self energy of the system in the 
presence of the measuring field. We first check that our approach reproduces correctly known results 
in simple limits, like the FCS of the resonant level system (AIM without Coulomb interaction). We 
then proceed to study the FCS for the AIM both perturbatively in the Coulomb interaction and in the 
Kondo regime at the Toulouse point (we also study a related model of a spinless single-site quantum 
dot coupled to two half-infinite metallic leads in the Luttinger liquid phase at a special interaction 
strength). At zero temperature the FCS turns out to be binomial for small voltages. For the 
generic case of arbitrary energy scales the FCS is shown to be captured very well by generalisations 
of the Levitov-Lesovik type formula. Surprisingly, the FCS for the AIM indicates a presence of 
coherent electron pair tunnelling in addition to conventional single-particle processes. By means of 
perturbative expansions around the Toulouse point we succeeded in showing the universality of the 
binomial FCS at zero temperature in linear response. Based on our general formula for the FCS we 
then argue for a more general binomial theorem stating that the linear response zero-temperature 
FCS for any interacting single-channel set-up is always binomial. 

PACS numbers: 72.10.Fk, 71.10.Pm, 73.63.-b 



I. INTRODUCTION 

The Anderson impurity model is one of the best stud- 
ied models in condensed matter theory 0, Q ■ Despite 
being exactly solvable by means of the Bethe Ansatz 
(BA) method in the wide range of equilibrium param- 
eters 0, 0, , its non-equilibrium properties are not yet 
fully understood. Notable exceptions are the works on 
the non-linear I — V characteristics 0, S Hi- It has 
first been realised by Schottky 0, that the current au- 
tocorrelation spectrum (sometimes also called noise) car- 
ries information about the charge of particles participat- 
ing in transport. The inve stig ation of these properties 
has been started recently |lfi ITTL IT^ . However, the 
current- voltage characteristics and noise spectra are only 
the lowest order moments of the full current distribution 
function, which is needed to completely characterise the 
transport properties of the system. Although it is still 
quite challenging to access even the noise correlations in 
experiments, in recent years it became possible to mea- 
sure the third irreducible moment (third cumulant) of 
the current distribution function |l3j . It turned out to 
carry information about the influence of the electromag- 
netic environment on the transport through the system 
under consideration [l4|. Moreover, it has was been ar- 
gued that the third cumulant is more suited for measur- 
ing the charge of current carrying excitations than the 
noise correlations [l5j . In order to meet future experi- 
mental needs it is therefore natural to analyse the full 
counting statistics of the AIM. 

The principal question is this: what are the effects of 



the electron-electron interactions on the FCS? Is it pos- 
sible to gain insight into the properties of a strongly cor- 
related electron system by studying its FCS distribution 
function? We provide at least a partial answer in this 
paper. The answer turns out to be on the negative side, 
though it is a constructive one: we find that the bino- 
mial statistics is universal in the low temperature linear 
response limit. (Interactions only affect the magnitude 
of the effective transmission coefficient.) At high voltage 
(temperature) the effects of the interactions are indeed 
profound (see main text) if more model dependent. 

The AIM model is characterised by a number of dif- 
ferent parameters: the electronic tunnelling amplitude 7 
between of the impurity level (which we also shall some- 
times call 'dot' later) and the external electrodes, its en- 
ergy Ao , and the strength of the Coulomb interaction on 
the dot U . There are three different transport regimes: 
(i) resonant level case, when U is vanishingly small in 
comparison to all other energy scales; (ii) Kondo dot 
regime, when U is large and when the dot level lies deep 
below the Fermi energies in the electrodes; (iii) mixed 
valence regime, which comprises all other possibilities. 
The most interesting situation is (ii) when the dot is 
permanently populated by a single electron. The so- 
called Kondo-resonance (also known as Abrikosov-Suhl 
resonance) in the local density of states leads to a signif- 
icant increase of conductivity, which has recently been 
observed in experiments on ultra-small quantum dots 
[Til Il7j . This phenomenon is a signature of the Kondo 
effect and is a result of the exchange interaction between 
the local spin degree of freedom on the dot with those 
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in the leads. One important feature of this effect is the 
fact that it is growing stronger as the temperature is low- 
ered. From the mathematical point of view that moans 
that the exchange interaction is a relevant operator in 
the renormalisation group (RG) sense, resulting in a new 
ground state where the local spin is absorbed, the leads 
are coherent and the conductance is maximal (perfect). 

The equilibrium Kondo model, being one special case 
of AIM model, is integrable by means of the BA tech- 
nique 0,^3- While it is possible to infer the non-linear 
I — V from the knowledge of scattering matrix, it is not 
yet clear whether an extraction of noise or any other 
higher order correlations of current is feasible. It has 
been pointed out by Toulouse > that the Kondo model 
allows a trivial diagonalisation for one special parame- 
ter constellation, when the whole Hamiltonian becomes 
quadratic in fermion fields. In addition to this very use- 
ful feature the Kondo model at the Toulouse point turns 
out to be representative for the low energy behaviour of 
a generic Kondo model [2(| , reproducing all essential de- 
tails of the latter in the low energy sector. For the Kondo 
dot a similar procedure has been developed by Emery and 
Kivelson in [2l| and refined by Schiller and Hershfield @ 
in order to access the non-equilibrium current-voltage as 
well as noise properties. As has been shown in p3 |. this 
approach can be applied to access the FCS as well. 

Contrary to the single-channel Kondo model, which 
maps on the conventional non-interacting resonant level 
(RL) model at the Toulouse point, the Kondo dot under 
the same conditions is described by the Majorana RL 
Hamiltonian [2(| ■ It has been demonstrated in [23L l24l| , 
that an RL coupled to two half-infinite Luttingcr liquids 
(LL) at the special interaction parameter g = 1/2 can 
be re-written in terms of a Majorana RL model as well. 
In this way one can obtain the exact FCS for a genuine 
interacting svstem. Related systems have been analysed 
before, see [2^,|2^| ('Coulomb blockade' dots). Below we 
investigate how the two models are related. 

The outline of the paper is as follows. In Section ITT1 we 
present a further development of the Levitov-Rcznikov 
(15l | approach to the FCS calculation in tunnelling set- 
ups. In order to test our Hamiltonian formalism we per- 
form an explicit calculation of the generating function 
for the FCS of a simple tunnelling contact between two 
metallic electrodes, see Section III Bl Section IIIII starts 
with a reproduction of the known FCS for the RL model, 
which is the spin-less version the of AIM set-up without 
Coulomb interaction. Next we derive Eq. JSUJ), which is 
the general formula for the FCS of an interacting sys- 
tem and the main result of this paper. We then proceed 
to evaluate the perturbative corrections in U 7 see Sec- 
tion IIII Dl and investigate linear response FCS on gen- 
eral grounds in Section IIII El The opposite case of large 
U, when the system is in the Kondo regime, is the subject 
of Section IIII Fl We not only present analytical results 
at the Toulouse point, but also analyse the change in the 
statistics around it in Section Till F 21 Subsequently we 
discuss the relation of the Kondo FCS to that of a RL 



set-up between two LL at g = 1/2 and establish connec- 
tions to existing results. Some conclusions are offered 
in Section IIVI There are several appendices containing 
technical details of some of the lengthier derivations. 



II. KELDYSH METHOD FOR THE 
CALCULATION OF CURRENT STATISTICS 

A. General considerations 

The cumulants of a given distribution function are 
known to define the latter in the unique way |27| . For 
practical reasons it is usually more convenient to calcu- 
late the so-called generating function x(A), which in case 
of charge transport is given by x(A) = J2 q etqX Pqi where 
P q is the probability for the charge q to be transferred 
through the system during the measuring time T. The 
parameter A here is referred to as the measuring field. 
The cumulants (S n q) (which are nothing else but the ir- 
reducible moments of P q ) can then be found for according 
to the prescription 

gn 



(S n q) = H) n Irxx(A) 



A=0 



The measurement of the charge transmitted through a 
system is usually accomplished by a coupling to a 'mea- 
suring device'. In the original work by Levitov and Leso- 
vik it is a fictitious spin- 1/2 galvanometer coupled to the 
current |2^, [2|j . The transmitted charge is then propor- 
tional to the net change of the spin phase. As has been 
shown by Nazarov [30| , the counting of charge can in gen- 
eral be done by coupling the system to a fictitious field 
and calculating the non-linear response, which leads, of 
course, to exactly the same results. 

According to |l5j the generating function is given by 
the following average, 



X(A) = ( T c exp 



-» / T x (t)dt 



(1) 



where C is the Keldysh contour, Tc is the contour order- 
ing operator, X(t) is the measuring field which is non-zero 
only during the measuring time T: X(t) = \9{t)6(T — t) 
on the forward path and X(t) = —\9(t)6(T — t) on the 
backward path. Introducing the operator transferring an 
electron through the system in the positive direction (i. c. 
in the direction of the current) Tr, and its counterpart 
Tl we can write 



JX(t)/2 T 



-i\(t)/2rp 



L ■ 



(2) 



We note in passing that T^ R = Tl in any system. Conse- 
quently, writing out Q explicitly in terms of the time- 
ordered and anti-time-ordered products, one arrives at 
the conjugation property 



X*(A)=X(-A) 



(3) 
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We now allow X(t) to be an arbitrary function on the 
Keldysh contour, X T (t) on the forward/backward path. 
Then a generalised counterpart of Q can be defined as 



(4) 



Next we assume that the measuring field changes only 
very slowly in time. Then up to the switching terms 
(which are known to be proportional to In T) 



x[A_(f),A+(f)] =cxp 



T 

i I U[X-(t),X+{t)]dt 
Jo 



where U(X-, A + ) is the adiabatic potential. Once the adi- 
abatic potential is computed, the statistics is recovered 
from 

lnx(A) = -iTU(X,-X) . 

Alternatively we can level off the X± functions in Eq. (@J 
to different constants as 

x[A_(i),A + (t)]->x(A- ) A+) , 

then x(A) = x(A, ~ A). Note that the conjugation prop- 
erty J2J) now generalises to 

X*(A_,A + ) = X (A+,A_), 

or 

W*(A_,A+) = -W(A+,A_) . 

To calculate the adiabatic potential we observe that ac- 
cording to the non-equilibrium version of the Feynman- 
Hcllmann theorem HH Ha, Hal , 



d 
dX~ 



W(A_,A+) 



dTxjt) \ 
9A_ / \ 



(5) 



where we use notation 



-ijT x (t)dt 

T G {A(t)e o 



X(A-,A+) 



(and similarly for multi-point averages) where A's are 
understood to be different constants on the two time 
branches. Note that the above one-point averages de- 
pend on the branch the time t is on (though not on the 
value of t on that branch): 

Therefore the average in Eq. JSJ) must be taken on the 
forward branch of the Keldysh contour. One immediate 
advantage of our Hamiltonian approach is the fact that 
the calculation of the adiabatic potential U amounts to a 
calculation of some well defined Green's function (GF), 
even though a non-equilibrium one. So we can use the 
whole power of the diagram technique and connect to 
many known results within this method without being 
restricted to scattering formalism as in [TM l29l | . 



B. FCS of a tunnelling junction 

In order to illustrate the procedure we calculate the 
FCS of the tunnelling junction between two metallic elec- 
trodes, denoted by R and L, which we model by the 
wide flat band Hamiltonians Hq[iPr 7 i,]. Their chemi- 
cal potentials are assumed to be prx = ±V/2, where 
V is the voltage applied across the junction (we set 
e = m . = h = ks = 1 and the Fermi energy Ep = 
throughout). The coupling between the electrodes is sup- 
posed to be the conventional point-like tunnelling with 
the amplitude 7, so that (for simplicity we assume spin- 
less electrons) 



H= fi oM + 7^Wi(0)+h 

i=R,L 



The unperturbed GFs (for 7 = 0) can be easily evaluated, 
see e. g. [H (i = R, L), 



9i M = 9^ + (oj) = i2Trp [n l 
g^ + {w) = i2ir p rii , 



1/2] 



■i2np [l 



(G) 



where pa is the density of states in the electrodes in the 
vicinity of Ep . Here ur^l = np(uj±V/2) where np is the 
Fermi distribution function. We use the original notation 
of Keldysh for the GFs, where the superscripts stand 
for the position of the time arguments on the contour C 
rather than the far more widespread language in terms 
of retarded (advanced) and thermodynamic components 
(EH HH vs. HI HI). The reason for that is the fact 
that due to the presence of two different fields A± the 
fundamental relation connecting the four Keldysh GFs, 
+ G ++ = G ^ + G" 1 , does not hold any more. 
Therefore in the present situation there are indeed four 
independent GFs. 

For obvious reasons the T\ operator Q is given by [l5j 



T\ = 7 



^(O)V'l(O) + e- iA Vl(0)Vfl(0) 



so that we have to evaluate 
d 



<9A_ 



W(A_, A+) = i 7 (e lA - j>y L - e~ iX - ^r)x 



(7) 



Defining the mixed GFs, 
G RL (t,t') = 



-i{Tci) R (t)i>l{t'))x 
GLRit.t 1 ) = -i(T c iP L m R (t'))x 



we can re- write (Q) as 
d 



9A_ 



W(A_,A+) 



lim / ^ 

£ ^o+ ./ 2n 



x [7*e iA -G^( W )- 7 e- <A -G^(a;) 



(8) 
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The calculation of the GFs G LR ^ RL ) (^) is most elegantly 
accomplished using functional integration. To that end 
we introduce the matrix of GFs according to 



9rr 9rr 9rl 9rl 

9rr 9rr 9rl 9rl 

9lr 9lr 9ll 9ll 

9lr 9lr 9ll 9ll 



Using © one easily constructs the corresponding matrix 
go without tunnelling. The GFs for 7 7^ are found from 
the equation 



G- 1 = g^ 1 -t t . 



(9) 



where N = TV/2ir = TVe 2 /h is the number of incoming 
particles during the waiting time (also known as 'number 
of attempts') and 1— T(0) and T(0) are their probabilities 
to be reflected or transmitted, respectively. Generally, 
the terms proportional to (e jmA — 1) may be interpreted 
as describing the tunnelling processes of particles with 
the elementary charge me [28| . Negative m correspond 
then to transport in direction opposite to that of the 
applied voltage. Due to the detailed balance principle, 
such terms do not contribute at T = 0. 



III. FCS OF THE ANDERSON IMPURITY 
PROBLEM 



where is the self-energy due to tunnelling. It is found 
to have only four non-zero components, 



t^rl 

T^LR 



7e 



iX. 



V+ + 
T RL 



iA-i 



7 e t^ lr 



-7 e 



Solving 10 results in 



Det (g- 1 



„u_ 7 



(^Po) 2 



[1 + T + 2(n R - n L )] , 
Det(^ - T t)G- LR =-e^-^ 



x [l + r + 2(l-2e iA )(n fl -n L )], 

where A = A_ — A + and T = (irpoj) 2 is the dimcnsionlcss 
contact transparency. The determinant is found to be 
given by 

Det (g^ 1 - T t) = (^o)" 4 {(l + T) 2 + 4r 
x [(e* A - l)n L (l - n R ) + {e,- fx - l)n R (l - n L )] } • 

Inserting this outcome into JSJ, integrating over A_ and 
setting A_ = — A + = A we immediately arrive at the 
Lcvitov-Lesovik formula 



ln Xo (A;V;{T(a,)}) = T 



2^ 



InU 



{l + T(w) (10) 
-1)]}- 



x [n L {l - n R )(e lX - 1) + n R (l - n L ){e 



where the transmission coefficient is given by T(lu) = 
4r/(l + r) 2 for the particular case of the tunnelling 
junction set-up. Eq. ijTU|) holds of course, for any non- 
interacting system, with known transmission coefficient, 
coupled to two non-interacting reservoirs described by 
filling factors n R ^. 

The generating function (jTU|) leads at zero temperature 
and small voltage to the conventional binomial distribu- 
tion function 



x(\)= [i-r(o) + r(o)« 



N 



A. Preliminaries 

Now we are in a position to proceed to more compli- 
cated models. The Hamiltonian of the AIM model con- 
sists of three contributions, 

H = Hq + Ht + H c ■ 
The kinetic part 

H = J2 H o[^R/L,a] + X)( A o + oh)dld a , 

describes a single fermionic level (which we shall also call 
'dot') with electron creation operators d\ (a is the spin 
index), energy Ao and subject to a local magnetic field h. 
Two non-interacting metallic leads i — R, L are modelled 
as in the previous Section. The leads and the dot are 
coupled via tunnelling, 



Ht = Y. [7Le 4A(t)/2 4^ 



iR^Rada^ 



H.c. 



with different amplitudes Tax. For convenience we al- 
ready included the counting field into the Hamiltonian. 
Notice that since the transfer of a physical electron 
through the device is a two-stage process (left lead — > 
dot — > right lead or the other way round) the measuring 
field is halved. For the sake of simplicity we incorpo- 
rate the counting field only into the left junction. Do- 
ing that at both junctions (of course with the correction 
A/2 — » A/4) leads to exactly the same results due to the 
gauge symmetry of the Hamiltonian. Finally, we include 
the Coulomb repulsion on the dot, 

H c = Un^n l , 

where n a = d\d a . The applied voltage is incorporated 
into the full Hamiltonian as in the previous Section, [Al — 
Hr = V> 0. 

We start with the definition of two auxiliary GFs, 

F x (t,t') = -i(Tc{Mt)dHt')}) x 
F x (t,t') = -i(T c {d(t)4,tf)}h 
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Hence the derivative of the adiabatic potential is given 
by 



d 



■U(X-,X+) 



1L lim / 



o+ / 2tt 



(11) 



Similar to the situation of the tunnelling junction these 
mixed GFs can be written as combinations of bare lead 
GFs and exact impurity GF D(t,t), 



F x (t,t') = / dt"g L (t-t") 



«\ -*A(t") nf .H .1 



D(t",t'), 



F x (t,t') = J dt"D(t,t')e- iX ^g L (t" -t') . 
c 

Performing the Kcldysh disentanglement and plugging 
the result back into one obtains 



U(X-,X+) (12) 



d 
dXl 

jl f dw [■ 
2tt 



where again A = A_ — A + . Hence, the whole problem 
is now reduced to calculation of the impurity GF. The 
most compact way to access it is using the self-energy 



J 



formalism. According to [3|j the self-energy E(w) for a 
non-equilibrium system can be defined in very much the 
same way as in the traditional diagram technique via 



D(u) = d (oj) + D(cj)E(u>)d (u>) , 
where the unperturbed dot GF is 



4 1 = 



w - A 
-w + A 



(13) 



(14) 



Then trivially 

b- 1 ^) = d^iuS) - ±{uj) . 

Therefore our goal now is the evaluation of the self- 
energy. 

B. The U — case: resonant level model 

We shall elaborate on the formula Ijl2|l . which is still 
valid for the interacting case, in the following subsection. 
For pedagogical reasons we pause here to deal with U = 
case, when H is trivially diagonalisable. This situation 
is referred to as the resonant level (RL) model. 

The corresponding self-energy is (we neglect the spin 
index here as GFs are diagonal in a and independent of 
it, the sub-script '0' distinguishes the [7 = quantities): 



E (u>) = 



1 2 l9l + 1r9r 

_ e -iX/2 7 2 g +- _ 7 2 5 +- 



iT L {2n L - 1) + iT R (2n R - 1) -2ie f ^ 2 T L n L - 2iT R n R 

2ie- lX / 2 T L {\ - n L ) + 2iT R {l - n R ) iT L {2n L - 1) + iT R {2n R - 1) 



where, in order to unburden the notation, we set T R j, = {it po') r ,l) 2 ■ Consequently 
D "V) = 



u - A - iT L {2n L - 1) - iT R (2n R - 1) 
-2ie- lX / 2 T L (l - n L ) - 2iT R {\ - n R ) 



2ie lX / 2 T L n L + 2iT R n R 
uj + A - iT L (2n L - 1) - iV R {2n R - 1) 



Inversion of this results in 



©o(w) 

w - A + iT L (2n L - 1) + iT R (2n R - 1) 
-2ie-^/ 2 r L (l - n L ) - 2iT R {\ - n R ) 



2ie lX / 2 T L n L + 2iT R n R 
-u + A + iT L (2n L - 1) + tT R (2n R - 1) 



where (r = r H + r i ) 

V (u) = (oj- A ) 2 + T 2 + 4r L r fl [n L (l - n R )(e* x / 2 - 1) + n R (l - n L ){e~ fx ' 2 - 1) 
Inserting these results back into (|12f) yields an equation, 



d 
<9A_ 



■U{X-,X+) = -2T L Y R 



duj n L (l - n R )e lX / 2 - n R {\ - n L )e~^/ 2 
2tt T> (w) 



(15) 



(16) 



(17) 



(18) 
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Performing the integration over A_ and constructing the generating function we again find the formula (|10|> with the 
Breit-Wigner transmission coefficient 



as expected for the RL set-ups |57j . 



T(w) 



w L r R 



(w- A )2 + r2 ' 



C. The general formula 

The GFs can be used to construct the consistent expansion of the FCS to all powers in U, opening the road to 
perturbative as well as non-perturbativc studies of the FCS. From now on under E we shall understand the self-energy 
due to the Coulomb interaction (tunnelling terms are incorporated into the bare GFs). Eq. 1)15(1 thus changes to 



ui - A _- iT L (2n L - 1) - iT R (2n R - 1) - E~ 2ie a / 2 Y L n L + 2iT R n R - £" 

-2ie-^/ 2 r L (l - n L ) - 2iT R (\ - n R ) - £+" -w + A - iT L (2n L - 1) - 4^(2/1^ — 1) — E H 



(19) 



After the inversion of this matrix and insertion it into (|12|> one gets [D{u>) is the corresponding counterpart to l|17(l ] 



d 



■M(A_,A+) = -r L 



{2r fl [e^ 2 n L (l - na) - e- 4X/2 n K (l - n L ) 



2-kV(uj) 

— OO 

e a / 2 n L V + ~ + e~ l ^ 2 (l - n L )IT+\ } 



(20) 



which is a general formula for the statistics in interacting systems. Here T>{u>) is the, A-dependent, determinant of 
the matrix given by Eq. (|19|) . For A = the rhs of this relation is proportional to the current through the device. 
Moreover, as expected, in this particular case Eq. I|2l)[l can be brought into the form derived by Meir-Wingreen [3^, 
when the transport is defined solely by the retarded dot level GF after a symmetrisation procedure. The presence of 
the counting field docs not allow a similar reduction for arbitrary A though. 

Clearly formula l|2(J|l is not restricted to the AIM as such but is applicable for any similar one-channel impurity 
set-up (including, e. g. electron-phonon interaction on the dot or a double dot). 



D. Perturbative expansion in the Coulomb interaction 



The obvious way to proceed is to calculate the lowest-order contributions to the self-energy, in the time domain 
E(t) = 



-iUD o -(0) + U 2 [D^-(t)] 2 D^-(-t) 
-U 2 [D+-(t)] 2 D + (-t) 



iUD+ + 



(t)] 2 D+~(-t) 



(0) 



U 2 [D++(t)} 2 D++(-t) 



The linear in U part is diagonal and is essentially a remnant of the occupation probability of the dot level (cftd). 
It is most conveniently evaluated in the following way (from now on we consider a symmetrically coupled system 
T R = Tl = r/2 at zero temperature in order to simplify the algebra) 



iUDr(o) 



-iUDo + 



(0) 



-iU 



/du> \ , 



where the object 



(1 + e iA ) 



tan 



A + V/2 



tan 



(A ± V/2)e 



-i\/2 



is, in general, A dependent. For Dq + see Eq. (tTB|) . Here n simply gives the dot occupation probability. Plugging 
this result into (|19f) and proceeding to l|20|l we find the result identical to 1|18[) up to the denominator l|17|) where the 
bare level energy Ao now gets renormalised, Ao — > Ao + Un\. Subsequent expansion in U and integration over energy 
results in a well controlled contribution which vanishes for the case of the symmetric Anderson model Ao = —U/2, 
to which case the following considerations are restricted. 
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We concentrate now on the correction at the second order in U. One way to access the self-energies is through the 
evaluation of the corresponding susceptibilities. We define them as [note the sub-script '0', not to confuse with the 
generating function x(^)] 



dui 
2^ 



Xo(fi) = i 

The respective self-energy can be extracted from 

S(w) = i 



D Q -(Lu + n)Dv-(uj) D^ + (lu + Q)D+-{lo) 
D+~{oj + tt)D + {cj) D+ + {uj + n)D++(u J ) 



2^ 



D+~(lu - Q)xo + m D+ + {lu - n)xt + m 



The equilibrium results have been originally presented in famous series of papers by Yosida-Yamada |40t l4ll |42| (we 
set T = for simplicity) , 



Eeq(w) = (1 - Xo)w 



1 

-1 



2T 



sign (a;) 26(—u>) 
—20(lu) sign(w) 



where the exact even-odd susceptibilities possess the following expansions in powers of U, 

.2\ V 2 jj 



Xc = 1 + ( 3 - — 



r 2 r 2 



+ 



Xo ~ Trr 



For finite V and at the second order in U, there are three distinct energy regions contributing to Eq. H20|) : —V/2 < 
u> < V/2, V/2 < lu < 3 V/2, and —3V/2 < u> < —V/2. The low-energy expansion (not only small U but small V as 
well) in presence of A one finds in the region — V/2 < lu < V/2: 



S(w) = (1 - Xo)w 
iU 2 



87r 2 r 3 



1 o 

o -1 

,-2iX ISV 



6uuV 



-i\ (3V 



-3e 



6luV 



Needless to say, these relations are consistent with the non-equilibrium calculation by Oguri j^j . On the other hand, 
for lu > V/2 one obtains 



-+(cu) = 



ie~ iX U 2 (W 
87r 2 r 3 



while for lu < V/2 the relation 



-.+- 



(lu) 



-2i\jj2 



U 2 fSV 



?7T 2 r 3 



holds. These self- energies, being incorporated into (|2(jp. yield the following generating function for the FCS, 

TTI 2 V 3 TTI 2 V 3 
lnx(A) = lnxo(A) + (e^ A - 1) + T7 ^ T ( e - 2 > A - 1) + 0(U A ) , 



where 



V/2 



i27r 3 r 4 



lnXo(A)=T J Ijln 



-V/2 



x 2 w 2 + r 2 



i) 



still contains U. Performing the expansion around the perfect transmission (hence the sign change of A in the following 
formulas) we see that in terms of susceptibilities 



lX + W 2 ~ 



Xe + Xo( e -iA _ ^ + 2^ e -2iA _ X j 



(21) 
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where N — TV/ir is the number of incoming particles during the measuring time slice. This is, of course, only valid 
at the order U 2 . We speculate that the general formula for the full FCS could be written in terms of the equilibrium 
susceptibilities Xo,e only. One possibility is the generating function of the form 



x(A)=ATln 



1 - 



i2r 2 



(e^-l) + ^V 2 (e^ x -l) 



(22) 



as this expression reproduces the expansion Eq. (|21|l . We stress again that so far we have only shown that Eq. 
holds at the second order in U and beyond that it is a mere hypothesis. 

It is tempting to interpret the appearance of the double exponential terms as an indication of a coherent tunnelling 
of electron pairs (caution: similar terms would also appear for the non-interacting RL model due to the energy 
dependence of the transmission coefficient). In the Toulouse limit calculation below we find further evidence for such 
interpretation. 



E. Linear response FCS 



rive at a conclusion that 



Here we would like to take a closer look onto the gen- 
eral formula (|20|l at zero temperature and vanishing ap- 
plied voltage. In order to arrive at correct results one has 
to bear in mind that the limits V — > and to — > do not 
commute in the presence of the counting field. Indeed, 
calculating the Keldysh determinant in both limits we 
see that 



lim lira V (u, V, A) = A 2 , + r 

uj—>Q V— >0 



(23) 



but 



lim lim VqUj.V.X) = A 



T 2 + AT L T R {e tX -l) . (24) 



In fact, it is the second scheme we have to implement 
analysing the first term in Eq. H2U|) . This leads to a trans- 
mission coefficient type contribution to the generating 
function. 

On the contrary, in the second term in Eq. I|2(J|) , which 
is produced by the self-energy, not even the integration 
over ui is restricted to [0, Vj . As a matter of fact, due 
to Auger type effects jij, |45j one expects that there 
are contributions to the current (and FCS) at all en- 
ergies. This effect is itself proportional to the applied 
voltage though, and results therefore in non-linear cor- 
rections to the FCS. Hence the energy integration can 
be regarded to be restricted to [0, V] even in the second 
term in Eq. 1|20|) . Moreover, since the self-energy does 
not have external lines and all the internal frequencies 
have to be integrated over, the limits V — > and ui — > 
in this case commute. That means that for the evalu- 
ation of the self-energy to the lowest order in V one is 
allowed to use the equilibrium GFs, calculated in pres- 
ence of the counting field A, i. e. (|16|l with tir = rij, = Uf 
and with the corresponding Keldysh denominator (|23fl . 
Therefore all diagonal Keldysh GFs are equal to those in 
the equilibrium and all off-diagonal ones are simply pro- 
portional to the same diagrams as in equilibrium. Since 
any given off-diagonal self-energy diagram describes an 
inelastic process, it should vanish for to — > and we ar- 



lim S(w) = Re£ rt (0) 

uj— ^0 



1 

-1 



even at finite A. Eq. (|20|) thus leads to the fundamental 
result 



lnx(A) = Nlnl 1 + 



r- 



(e a -lU ,(25) 



[Re£(«)(0)] 2 +r 2 



or to In x(A) = iXN for the symmetric Anderson impurity 
model. In case of the asymmetrically coupled impurity, 
T/f 7^ Tl the numerator of l|25[) modifies to r^r^ while 
the denominator contains (r# + Tl)/2 instead of T. 

The result <|25|) allows simple generalisations to asym- 
metric svstems in a magnetic field h. According to 
|4oL l4lL E^ j the real part of the self-energy is given by 

Re4 i?) (0) = XcK + a Xs /i, 

where x c / s are exact charge/spin susceptibilities (combi- 
nations of even/odd) and k ~ Ao + U/2 is a particle-hole 
symmetry breaking field. Consequently 



lnx(A) = — In 



1 



1 



[xcK + xsh] 2 + r 2 
r 2 



(e iX 1) 



(e lA - 1) 



.(26) 



The enormous advantage of this formula is the fact, that 
the susceptibilities can be calculated exactly for any sys- 
temparameters with the help of the Bcthe-Ansatz results 

Let us stress that the result (12511 is not limited to the 
AIM but will hold for any similar model, hence the bi- 
nomial theorem. It is clear in hindsight that all the non- 
elastic processes fall out in the T = linear response 
limit . Still it is a remarkable result that all moments have 
a simple expression in terms of a single number: the effec- 
tive transmission coefficient. The binomial distribution 
is universal. [For a multi-channel system modifications 
will be required as is obvious from looking at Eq. Q26[)]. 



J 
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F. The Kondo regime 

The way to proceed further is to consider the case of very deep Ao and strong Coulomb repulsion. In this limiting 
case the system is in the Kondo regime and the dot can in good approximation be considered to be permane ntly 
populated by a single electron. It has been shown in [4^, that the conventional Schrieffer-Wolf transformation |47j . 
which maps the Anderson impurity Hamiltonian onto that of the Kondo problem, also works out of equilibrium. The 
result is the two-channel Kondo Hamiltonian 

H = H + Hj + Hy + Hm , 

where, with tp a ,tr are the electron field operators in the a = R,L electrodes, 



Hn 



i i J2 E / dx^ aa (x)d x ^( X ) , Hj= J2 E 



a=R,L <t=1,1 



H V = (V/2)£ J dx&l^r 



a,/3=R,L v=x,y,z 

H M = -\±B9ihT z = -At 2 



(27) 



Hb is the Bohr's magneton, gi the gyromagnctic ratio and h denotes the local magnetic field, which is applied to the 
impurity spin. Here t u=x * v,z are the Pauli matrices for the impurity spin and 



^=E^-(°) ct -'^'(°)' 



are the components of the electron spin densities in (or across) the leads, biased by a finite voltage V . The last term in 
Eq. (|2Z|) stands for the magnetic field, A = f i B g l h. We follow and assume = = jf , J z ± = (J^ L ±J RR )/2 
and Jg R — J RL = 0. The only transport process then allowed is the spin-flip tunnelling (sometimes also called 
'exchange co-tunnelling'), so that we obtain for the T\ operator 



J 



RL 



-V^/^Vlt +T-e iA »/VV^ +r + e~ iX ^t L i> m 1^%^) 



Of course, there is also a regular clastic co-tunnelling term, which couples the leads directly. However, it can be 
rigorously shown |46[ , that these processes are subleading in the low energy sector in comparison to spin-flip tunnelling. 
That is why we keep only the latter contributions to the Hamiltonian. We proceed by bosonization, Emcry-Kivelson 
rotation, and refermionization [l,|2iJ,El. We obtain then with J± = (jf L ± J RR )/y/2lrao, J± = J RL /y/2-Ka (ao is 
the lattice constant of the underlying lattice model) 



H = i 



E 

I/=C,3,c/,s/ 

J 



dx^ v (x)d x Mx) +^ [^/(o) + ^/(o) 
'■Y - ^/(o)] + 3 -f [^(o) - v«/(o)_ 

( J z+ - 2vr) : ^1(0)^(0) : +J Z - : ^(0)^/(0) 



V / dxipl f (x)ip cf (x) . 



(28) 



where now four fermionic channels are present: (c) total charge density channel for the sum of particle densities 
in both electrodes, (cf) charge flavour channel for the difference in densities. The channel-symmetric spin density 
channel (s) and channel-antisymmetric (or spin flavour channel) (sf) (see details in || and po|') are defined in analogy 
to their charge counterparts. A considerable simplification of the theory is achieved by introduction of the Majorana 
components of the continuum fields 



(29) 



and of the impurity spin t x = b and t v = a. 
operators in the terms quartic in fermions) 



As a result, the model simplifies to (for convenience we kept ips,sf 

(30) 



H = H' + H" 
H' = H' -i(J-b£ sf + J + ar] s f)-iAab + T\ 



H h 



v 1 :^t(0)fl>,(0):+J,---il>lf (0)^/(0) 
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where v i = J z + — 2n and the counting term is given by 

T x = -U ± b %} cos(A/2) - ncf sin(A/2)] . 



(31) 



The fields r/ s f and £ s j in the spin-flavour sector are equilibrium Majorana fields, whereas rj c f and £ c f in the charge- 
flavour sector are biased by V (from now on we omit the cf index and denote sf by /), 



H' Q = i 



dx 



r)f{x)d x r) f (x) + £ f (x)d x £ f (x) + rj(x)d x rj(x) + £,(x)d x £,{x) +V£(x)r)(x) 



where we drop the c and s channels as they decouple from the impurity completely (at the Toulouse point). The 
evaluation of the adiabatic potential can now be performed along the lines of Section ITT1 



J-U(X-,\ + ) = -^ J g[sin(A_/2)G b yH+cos(A_/2)G- r H 

where we again define the mixed GFs according to the prescription 

Gtt(t, t') = -i{T c b(t)C(t))x , Gbr,(t, t') = -i(T c b(t) v (t)) x . 



(32) 



(33) 



1. The FCS at the Toulouse point 

For realistic systems it is reasonable to assume J z _ = 0. The only remaining term which is still quartic in fermionic 
fields is then zero at the so-called Toulouse point J z + = 2n [Tflll20j|. In this situation the Hamiltonian is quadratic in 
fermionic fields. The mixed GFs l|33|) are related to the exact impurity GFs, D^t^t') = —i(Tcb(t)b(t')) \ and to bare 
GFs (calculated for all Ji = 0) for the Majorana fields |24| (notice that in the present situation we have to double the 
applied voltage in comparison to ©, for details see |24|). 



3« - Qrm - o 



tir + til — 1 riR + til 
n.R + UL — 2 tlr + til - 1 



» 9n£ = 



n L - n R 



1 1 
1 1 



(34) 



in the following way, 

G K {t,t') = ^J^/dm^^t'OI-cosIA^O^a^-iO + sinlA^O/^^'-t')} , 

G bn (t, t') = U ± J dt"D bb (t, t") {- cos[A(t")/2].g^(t" - f) + sm[\(t")/2]g m (t" - t')} . 
c 

After the Keldysh disentanglement and using (13411 

gf-W(A_, A+) = «4 / ^ {£»»-M(n B - n L ) + £>-+M [e^ 2 (l - n R ) e^l - n L )] } 



(35) 



Evaluation of the impurity GF is accomplished by the calculation of the corresponding self-energy and inversion of 
the emerging matrix d^ 1 — Y^k where in the absence of the magnetic field A = and J + = (we discuss the general 
case later) 



'1% 



(36) 



where the super-script (0) distinguishes the equilibrium GFs for V = and c = cos[(A_ — A+)/2], s = sin[(A_ — A+)/2] 
and <i _1 is given in Ijl4(l with Aq = 0. Using Ij34(l and new definitions I\ = J 2 /2 we obtain 



-<K 



T^(2n F - 1) +T ± (n R + n L - 1) 



,-iX/2 



nR.) 



r_ 2(1 - n F ) + t_l 

Then the determinant 



e a /2(i _ nR ) + e -iA/2(i _ nL ) r_(2n F - 1) + T ± {n R + n L - 1) 



- Dot (do 1 - t K ) =lu 2 + (Tj. + T_) 2 + ri [n L (l - n B ){e lX - 1) + n R (l - n L ){e~ lX - 1) 
+ 2r_r ± {[n F (l - n R ) + n L {\ - n F )](e 1 ^ 2 - 1) + [n F (l - n L ) + n R (l - n F )]( e - ll/2 - 1)} 



(37) 
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The GFs of interest are then given by 
Dct (do 1 - ± K )Dbb = 

-u + i[r_(2n F - 1) + T ± (n R + n L - 1) 



r_ 2n F + T x {S/ 2 n L + e- lX ' 2 n R ) 



iX/2 



-i {r_ 2(1 - n F ) + r ± [e^ 2 (l - n R ) + e- fx / 2 {l - n L )]\ w + i[T_(2n F - 1) + T ± (n R + n L - 1)] 

Inserting the calculated GFs into the fundamental relation H35[) results in 

d 7/(\ \ or f dLU I( ^> 
— ZY(A T ) = -z2I\l J 

with 



2^ Det^o 1 -%k) 



(38) 



n L (l-n R )e l - n R (l - n L )e 



-iX 



2r_r 



n F (l- n R y~ x ' 2 -n F {l-n L )e- ix ' 2 



To proceed, we split the ^-integral in Eq. I|38|) into two parts for negative and positive energies and change lo 
the second integral. In doing so observe that under this transformation n F — > 1 — n F , n R — > 1 — ul, and nj, — 
Therefore the denominator stays invariant while the numerator changes as 



-> — lo m 
l-n R . 



I(-co) = Ti [7i L (l-n i? )e lX -n ii (l-n L )e^ 
Eq. thus becomes 



2r_r_L 



n L (l - n F )e lX / 2 - n fl (l - n F ) e - fx ' 2 



d 



■W(A, 



2 7 2tt Det^V 1 -t K ) 



Observe that, crucially, = %[I(lo) + I(— u>)] so that the A-integration can be performed as before. The following 
exact formula for the statistics, valid at finite temperatures, follows immediately: 



lnx(A) = T J g In {1 + T 2 (lo) [n L {\ - n R )(e 2iX - 1) + n R (l - n L )(e- 2iX - 1)] 



(39) 



+Ti(o;) [[n F (l - n R ) + n L {\ - n F )](e tX - 1) + [n F (l - n L ) + n R (l - n F )}(> 



where the effective "transmission coefficients" (two of them now) are: 

rj , > _ 2r_r x 

^ 2 + (r_ + r ± ) 2 ' 1(UJ> lu 2 + (t. + t ± ) 2 



-iX 



1)]} 



T 2 (u) 



(40) 



In the more general case of finite magnetic field A and T+ the result l|39() is exactly the same up to the modified 
transmission coefficients (derived in Appendix A), 

ri^ + r 2 ,) 



T, 



- r+(r ± + r_)] 2 + lo 2 {t + + r_ + r ± ) 2 
2r ± r_(cj 2 + r 2 + ) + 2A 2 r ± r + 

- r+(r ± + r_)] 2 + ^ 2 (r+ + r_ + r ± y 



(41) 



In fact, since the refermionised Hamiltonian describes local scattering of non-interacting (Majorana) particles, the 
result H39f) can as well be derived using the approach originally conceived by Levitov and Lesovik for systems with 
known scattering matrix |28j |. For the corresponding calculation see Appendix B. 

Using the properties n F (l — n R ) + ulO- — n F) = — n R )(\ + cxp[— V/T]) and n F (\ — Ul) + n R (l — n F ) = 
n R (l — Jii)(l + cxp[V/T]) we can rewrite the result in the form 



ln X (A) = T / ^ ln{l + n L {\ - n R ) [r 2 (uj)(e 2 * x - 1) + T 1 (c)(e jA - 1)(1 + e" v ' T ) 
o 

+n R (l - n L ) [T 2 (cu)(e- 2lX - 1) + T^lo)^ - 1)(1 + e v / T )] } , 



(42) 



We first take a look onto the T = situation, when cxp(— |V|/T) — ► 0. In that case one can reduce the gen- 
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crating function to the Levitov-Lesovik formula (|10[l for 
a spinful system |^ 



X (A)= [l + T e (e iX -l)] 



(43) 



where T e = y/T 2 (0) = r ± /(T ± + T_). Hence in the 
low temperature limit we obtain the conventional bino- 
mial statistics for the charge transfer through the dot. 
Needless to say this is in accordance with the binomial 
theorem stated in the previous Section. However, the re- 
duction H43JI is not possible for finite temperatures and 
voltages. To the best of our knowledge Eq. I|39|l is the 
first exact result showing non-trivial statistics at finite 
energy scales. It can be interpreted in terms of two dis- 
tinct tunnelling processes: (i) tunnelling of single elec- 
trons and (ii) tunnelling of electron pairs with opposite 
spins. As has already been realised in at least in the 



regime T, V <C A tunnelling of single electrons is ener- 
getically very costly as it requires a spin-flip. A simulta- 
neous tunnelling of two electrons, which is described by 
the terms with 2 A and T 2 (cj), leaves the dot spin effec- 
tively untouched, making that kind of process the domi- 
nant transport channel. In zero field the finite voltage is 
known to act as effective magnetic field §| so that this 
tunnelling mechanism is always present regardless of the 
precise value of A. 

In the low energy sector to <C V, T the integral of H42|) 
can be performed explicitly in the spirit of Ref . j2|j , re- 
sulting in 



X(X) = exp 



where v = V/T and 



cosh(w) — cosh(v) = Ti[cos A — 1 + cosh(t> + iX) — coshu] + T 2 [cosh(u + i2X) — coshv] . 
In the limiting case V>Twc recover the result (|43H while for VcTwc obtain 

X(A) = exp (-t|-A 2 

where sin (A*/2) = 4T e sin 2 (A/2) [l - T e sin 2 (A/2)] . The full transport coefficient T as calculated in turns out to 
be a composite one and it is recovered from Ti,2 through a very simple relation: To = T 2 + T\/2. We have evaluated 
the first and the second cumulant of the Kondo FCS Eq. (|39|l which are the same as calculated by SH at all V and 
T We shall not reproduce these two cumulants here and concentrate instead on new results. 

First we would like to analyse the equilibrium statistics at V = 0. From Q39|l it is obvious that as Ul = tir all odd 
order cumulants are identically zero. Then for the even order cumulants wc obtain 



<^ 2 } 


= T 


(Sq 4 ) 


= T 


(*3 6 ) 


= T 



/ 7T [4(Ti + 8T 2 )n F (l - n F ) - 48(Ti + 2T 2 ) 2 n 2 F {\ - n F ) 2 } , 
Jo 111 

/ 7T 4 ( T i + 32T 2 )n F (l - n F ) + 1920(Tj + 2T 2 fn F {\ - n F ) 3 
Jo 2tt L 

240(T! + 2T 2 )(T! + 8T 2 )n|(l - n F ) 2 ^ 



As for finite /(0) one obtains / °° dujf{w)n F {l - n F ) n « a„T/(0) with oi = 1/2, a 2 = 1/12, a 3 = 1/60, a 4 = 1/280 
etc. all equilibrium cumulants are linear in temperature in the low energy sector. The lowest order cumulant is then 
the conventional thermal Johnson-Nyquist noise Sjn ~ 4GoToT, where Go is the conductance quantum and To is the 
transmission coefficient of the dot at ui = 0. 

In the opposite limit of finite voltage and T = 0wc obtain for the third cumulant 

(5q 3 ) = T — [T x + 8T 2 - 3(Tj + 2T 2 )(Ti + 4T 2 ) + 2(Ti + 2T 2 ) 3 ] . 

Jo 27T 

This simplifies further in zero field: 

T f 9W 2 

(Sq 3 ) = — ^2T ± t a n- l [v/(r ± + r_)] - ^ ^ +v2]2 [(r± + r_) 2 + 2r_(r ± + r_) + w 2 \ 



i 

possessing the following limiting forms: At low voltages the cumulant is negative for T_ < r j 

(8q 3 ) y^^TGo 2r y; + rr)3 r±) V ' (44) 
(Sq 3 ) v^oo w T7rG r ± . 
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2 4 6 8 10 



FIG. 1: Zero temperature voltage dependence of the third 
cumulant for different magnetic field values and T±/F± — 0.1. 
Inset: temperature evolution of the curve for A/T± = 1.5 for 
T/T x = 0, 0.2 and 1.5 (from bottom to top). 

Generally, under these conditions the n-th cumulant ap- 
pears to possess n~ 2 zeroes as a function of V, according 
to numerics. The saturation value in the limit V — ► oo is 
independent of the coupling in the spin-flavour channel 
because the fluctuations in the biased conducting charge- 
flavour channel are much more pronounced than those 
in the spin-flavour channel, which experiences only rela- 
tively weak equilibrium fluctuations. 

For the general situation of arbitrary parameters, the 
cumulants can be calculated numerically. The asymp- 
totic value of the third cumulant at high voltages, simi- 
larly to the findings of j2|| , does not depend on tempera- 
ture and is given by the result (|44|l . see Fig. 1 of [22|. In 
the opposite limit of small V, (Sq 3 ) can be negative. Suf- 
ficiently large coupling T_ or magnetic field, see Fig. ^ 
suppress this effect though. 



According to the result of Ref. [ll| , as long as the dis- 
tribution is binomial, (Sq 3 ) / (Sq) = (e*) 2 , where e* is the 
effective charge of the current carriers. This quantity is 
to be preferred to the Schottky formula because of its 
weak temperature dependence. Indeed we find numeri- 
cally that the ratio (Sq 3 )/(5q) in the present problem is 
weakly temperature dependent (it is flat and levels off to 
1) in comparison to (5q 2 ) / (Sq) . 



2. Corrections around the Toulouse point 

Thus far we dealt with a system which finds itself at 
one special point in the parameter space, when v\ = 
and J z _ = 0. While the latter requirement is reason- 
able for realistic systems, the former is quite artificial. It 
has been shown by means of RG transformation proce- 
dure that at least in equilibrium the operators, describing 
deviations from the Toulouse point are irrelevant in the 
RG sense and do not influence the physics in the low en- 
ergy sector strongly. There is, however, no a priori rea- 
son why that should hold in a non-equilibrium situation. 
Therefore the full analysis of the FCS must incorporate 
the investigation of the statistics beyond the Toulouse re- 
strictions. We first concentrate on the situation of finite 
v\. As was pointed out above, an analytic solution in this 
situation is not possible. The only option to progress is 
perturbation theory in v\ . 

To access the generating function we still can use the 
fundamental relation l|35|) . As we have the complete 
knowledge of all GFs with respect to H', see <|3(J[) . the 
simplest thing we can do is to calculate perturbative cor- 
rections to D bb in the second order in V\ . They are given 
by 

Dbb = D bb + v\D bb £ b D bb = D bb + SD bb . 
The correction to the adiabatic potential is then given by 



dU 



du> 
2^ 



{SD^-(n R - n L ) + SD M + [e lX / 2 (l - n R ) - e -' iX / 2 (l - n L )] } 



The self-energy matrix components are defined as 



2- 



de 



B2,(w- ei ) / ^G^( £2 )Gf( ei +£ 2 ), 



where G s (t, t') = —i(Tc4>s{t)i>\{t'))\ is the GF of the spin sector fermion which is free. Therefore at T = it is easily 
found to be 



In the A = and T + = case only D aa < ' ++ ' = ±^o = ±l/w are non-zero, so that only S fe needs to be calculated, 
resulting in 



-|sgn(tj) i6(— uj) 
-i@(uj) — |sgn(w) 



E t — (w) = -E++M = (2n)- 2 co (In | w | - 1) 
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It is an odd function of u> and vanishes in the infrared limit as expected from RG arguments since it is generated by 
an irrelevant operator. The corrections to the impurity GFs are then 



SD M - =vlZ-- (D-D" -D-+D+-) , SD- + =vl^--D; b + (D- b --D^ 



Furthermore, 



n—n— n-+n+- - Uot^f^ v \ i + » F-(2n F - l) + r ± (n fl + » L - l)]} 

u bb v bb - u bb D bb - Uct \ d a - h i<) + ~ a ,3-1 — a - ; ' 

Det (a — uk) 

and is an even function of to. After multiplication with the self-energy, which is an odd function, and after an 
integration over all frequencies we immediately see, that the off-Toulouse correction to the first contribution in 135(1 
is identically zero: 

D bb (Dbb ~D bb ) = -2*0, Dct2 (V _ ±K) ■ 

Therefore for the correction to the derivative of the adiabatic potential we obtain 

3& - I s 5?fa) " < 2r -"' [e,s,2(1 - "*> - e_,s/2(1 - " L)i 

+T ± [e i V(l - n R ) - e- fx n R (l - 7i L )]} , 

Comparing this result with ((38(1 we conclude that the effect of V\ is the correction to the transmission coefficients 
Ti = Ti + 5Ti with 



ST, 



since schematically the structure of the correction is 

2T 2 e 2lX + T x e iX . 2T 2 e 2lX + T x e iX 



1 + T 2 (e 2iX - 1) + Ti(e iX - 1) [1 + T 2 (e 2 * x - 1) + Ti{e iX - l)] 2 ' 
which can be seen as the lowest order expansion of 

2r 2 (l + <5)e 2 ' A +r 1 (l + £)e' A 
1 + T 2 (l + 5)(e 2iX - 1) + Tj(l + 5)(e lX - 1) ' 

This correction vanishes for uj — *• 0, hence the trivialisation ((43(1 still holds for the transmission coefficients away from 
the Toulouse point. 

It is not clear whether this picture is valid for T + , A ^ 0. As in this situation the magnetic field couples the a and 
b Majorana fields, the D aa correlation functions have to be calculated from the Dyson equation with respect to the 

self-energy = (ij) A 2 -D!1(A = 0) + T+ij^, . As our ultimate goal is the behaviour of the generating function 

in the limiting case of small V, we set 1^ = 0. We return to the finite V situation at the end of this Section. The 
resulting can be best given in terms of the retarded component , 

Ep (++) = d=ReSf -i(2njr-l)ImSf , 

= i2n F ImEf , (45) 

— „' 0/1 ~ A T™\ 

J b ' 



= -«2(l-nF)ImSf 



with 



Ef = 1 ^ [O 2 + A 2 - (T_ + T ± ) 2 ] {(flj - iw) [ln(^- - iw) - 1] - fi^lnfi,- - 1)} 
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where 6= (T_ + T ± ) 2 + T 2 + - 2A 2 , d= [A 2 + r+(T_ + T ± )] 2 , and Vt\ 2 = (6± V& 2 - 4d)/2. The expansion for small 
energies is different from that at A = 0, 

Re £« a " V p 2 + A 2 - (T_ + T ± ) 2 } lnfy 



(2tt) 2 V& 2 - 4d .fr'. 



J 

2 



i mE * „ _ f " 2 y [n 2 + a 2 - (r_ + r ± ) 2 ]^ 



The correction to the time ordered part is, 



SD^- = v\ (D^-D^- + D^+ + D+- + Di b + J:+-D^- + D^+D+- 
= v 2 (Da-Da- - D;+D+-) RcSf + l v 2 lm^T 1 (u J ) , 



with 



= -(2n F - 1) {D--D-- + D-+D+-) - 2(1 - n^D^D^ + 2n P D„-D 



bb 



This function is an odd function of ui while D bb D bb —D bb + D bb is even. Taking into account the symmetry properties 
of the self-energy we conclude that the whole contribution to the generating function stemming from 5D bb ~ vanishes. 
The other component can be written down in the similar way, 

5D-+ = v 2 (^r^r V + + v^r V + o fc 7^ + 



WPtt" - D b b + )^b + **? ImSf T 2 {uj) . 



where we have introduced 



^aM = (1 - 2n F )D bb +(D bb - + D++) + 2(1 - n F )D bb + D bb + - 2n F D bb ~Dt b + . (46) 

The analysis of the contribution arising from ReS^ can be done in the same way as before as it has exactly the 
same structure. The substitution l|50|l still can be applied and one immediately recognises that it only leads to the 
renormalisation of the transmission coefficients, 

ST t (uj) = « 2 ReSf u^ + Tl-tf) 

[to 2 - A 2 - t+(t ± + r_)f + ^ 2 (r_ + r+ + r ± ) 2 

ReS^ is itself linear in to in the low energy sector, that is why the corrections to the transmission coefficients (|47|) 
vanish at low energies. Now we turn to the contribution of ImS^. The first term in (|46J) can be shown to produce 
renormalisation of the transmission coefficient similar to (|47|l . 

oTAuj) = WilmSf 5 Tdbj) . 

[uj 2 - A 2 - r + (r^ + r_)] + w 2 (r_ + r + + r^) 2 

Although the two remaining terms of 1)46(1 cannot be reduced to renormalisation of the transmission coefficients in a 
simple way, their contribution to the derivative of the adiabatic potential can be evaluated directly, 

* (IB - <4 [ t ' <*' - *«> [(» - £kf + (- + ^/ + ^ 

Taking into account that the leading behaviour of the imaginary part of the self-energy is ~ uj 2 for small energies 
one immediately verifies, that the above correction is cubic in the applied voltage and therefore leads to qualitatively 
the same picture as the renormalisation of the transmission coefficients. For that particular evaluation we used 
the equilibrium self-energy £&. Nevertheless, after a lengthy but straightforward calculation, we find that the same 
conclusion is still valid for the proper non-equilibrium one as the corresponding corrections to Ef, is of exactly the 
same order in V and uj. Thus at least in the low energy sector the predictions of Sec. IIII F II remain valid beyond the 
Toulouse point. 
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3. Resonant level problem in Luttinger liquids 



We now briefly turn to the g = 1/2 RL set-up. This set-up has caused much interest recently, see Ref. (23| and 
references therein. The Hamiltonian now is 

H = H + ( 7L i> L S + 1R d^ R + H.c) + AdU + H c , 

where Hq stands for two biased Luttingcr liquids (LLs), d is the electron operator on the dot, 7 .r(l) are the tunnelling 
amplitudes to R(L) electrode and He is an electrostatic interaction (see also p^L 

H c = X c Sd 5^(0)^(0). 

i 

The contacting electrodes are supposed to be one-dimensional half-infinite electron systems. We model them by chiral 
fermions living in an infinite system: the negative half-axis then describes the particles moving towards the boundary, 
while the positive half-axis carries electrons moving away from the end of the system. In the bosonic representation 
-f/olV'i] are diagonal even in presence of interactions (for a recent review see e. g. |2(i| ; we set the renormalised Fermi 
velocity v = vp/g = 1, the bare velocity being vp ): 



HobPi] = (4.)- J dx[d x ^(x)Y. 

Here the phase fields 4>i{x) describe the slow varying spatial component of the electron density (plasmons), 

i)\{x)ij)i{x) = d x <pi{x)/2-K^/g . 
The electron field operator at the boundary is given by j5(| , 

^(0) = e^(°)/v^/V2^, 

where ao is the lattice constant of the underlying lattice model. Here g is the conventional LL parameter (coupling 
constant) connected to the bare interaction strength U via g = (1 + Uj-Kvpy 1 / 2 [US- In the chiral formulation 
the bias voltage amounts to a difference in the densities of the incoming particles in both channels far away from 
the constriction [4flL IHfl | . The current is then proportional to the difference between the densities of incoming and 
outgoing particles within each channel. 

Construction of the operator © is unproblematic and leads to 

T X = 7L(dW V4 + ^Ide- 1 ^) + lR (d^ R e-^ + ^ R de iX / 4 ) , 

where, contrary to the Anderson impurity calculation, we choose to build in the counting field in a symmetric manner 
for the reasons which will become clear later. After the Emery-Kivelson rotation, refermionization to new fermions 
■0 and after the introduction of the Majorana components as in 1|29[1 [2~Lll24l | we find 

T x = [e lX '\ lL £yj + j R tpd) + e~ lA / 4 ( 7L dV + lR ^ d ) 

= -i 7+ b [cos(A/4)£ - sin(A/4)r/] + a [sin(A/4)£ - cos(A/4)t?] . (48) 

where 7 ± = 7 l ± 7 ^. In case of the symmetric coupling 7 _ = the corresponding T\ has exactly the same shape as 
In fact we find the same set of equations as for the Kondo dot, Eq. (|30|) and Eq. i|31|) . but with A — > A/2 and 
J± = 7 +, J± =0. Consequently, the FCS is given by a modification of the Levitov-Lesovik formula i|l(J[l : 

Xi/ 2 (A) =xo(\;2V;{T a (lo)}) , (49) 

with the effective transmission coefficient Ta(w) = 4 7 4 o> 2 /[4 7 4 ijj 2 + (lo 2 — A 2 ) 2 ] of the RL set-up in the symmetric 
case [23|. All the cumulants are thus obtainable from those of the non- interacting statistics Eq. 1|1(J|I . 

The A = RL set-up is equivalent to the model of direct tunnelling between two g = 2 LLs 24] . The latter 
model is connected by the strong to weak coupling (l/g — > g) duality arg ument to the g = 1/2 Kane and Fisher 
model [48ll5lT|. which is, in turn, equivalent to the CB set-up studied in [25|,[2(| (for a more general case of arbitrary 
interaction strength see also p2i l53|\ Therefore their FCS must be related to our Eq. H49(l at A = by means of 
the transformation: Tq — * 1 — To and V — * V/2. Indeed after some algebraic manipulation with Eq. (12) of j2||, for 
details see Appendix C, we find that the FCS for the CB set-up can be re-written as: 

XcbW =xo(-X;V;{1-T (lu)}) . 

For the asymmetric coupling 7 _ ^ the problem cannot be mapped onto the Kondo dot any more. The correspond- 
ing calculation is nevertheless straightforward and is presented in Appendix D. There is no fundamental difference 
in the result up to the more involved transmission coefficient, which has already been derived for the case of the 
non-linear / — V in |23l ■ 
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IV. CONCLUSIONS 



To conclude, we present a detailed study of the charge transfer statistics through the Anderson impurity model. We 
find an expression for the exact generating function in terms of the impurity self-energy calculated in the presence of 
the measuring field A: Eq. i|20[l . Based on this formula we conclude that T = linear response statistics is universal 
and binomial for the AIM and similar models: we call this fact binomial theorem. The only effect of correlations is to 
define an effective transmission coefficient. For the symmetric AIM, for example, there is a perfect transmission and 
no fluctuations of the current at all in this case. 

In the search for non-trivial interaction effects one has, therefore, to go to higher values of T and V. To this end we 
have calculated the exact FCS distribution function in the Toulouse limit (Kondo regime): it is given by Eq. I|39|l . This 
formula uncovers rather profound, if model dependent, consequences of correlations: there arc two distinct tunnelling 
processes (T± and T2), that of single electrons and electron pairs with opposite spin. The latter process is, in fact, 
dominant in zero field. The structure of higher moments is also determined by these two processes as discussed in 
detail in the main text. At T = linear response all this rich physics is masked by Eq. (|39|) collapsing to the universal 
binomial distribution. We checked this universality by extensively studying corrections to the distribution function 
due to departures from the Toulouse limit. 

We close by outlining some possible directions for future developments. Formula i|20|) could be used to develop 
Fermi liquid theory for the noise and possibly higher moments. Perhaps more importantly, the ideas of this paper 
could be applied to models with many conduction channels, where one would expect some equivalent of the binomial 
theorem to hold, as seems to be compatible with recent experiments |54j . 
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Appendix A 



The calculation of the impurity GF in finite magnetic field A and J+ 7^ is accomplished again by inversion of 
oIq 1 — Sa with 



E A = 



A 2 D: 



-A 2 D~+ + E-+ 
-A'£>+- + E+- A 2 £++ + E++ 



where is given in l|3tj|) . D aa has to be evaluated with respect to the Hamiltonian 

H+ = H [r)f] - iJ+arjf . 

A relatively simple calculation yields 



C0 2 +T 2 



uj + iT + {2np — 1) i2T + np 
-i2T-(l-n F ) -u + iY + {2n F - 1) 



From now on the calculation can be performed in exactly the same way as before. However, writing down explicitly 
the expression for £a one observes, that it can be constructed from the corresponding Sjf via trivial substitution 



A 2 u 
; 2 + T 2 + 



A 2 r, 



(50) 



This can be used to obtain the transmission coefficients (|41|) from the ones given by 1)40(1. 
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Appendix B 

For systems with known scattering matrix s a (3 imn between terminals a, (3 and channels m and n, there is a ready 
formula for the FCS generating function, derived in [2gj . 



Inx(A) = f- 



dio In Dot 



l + />) (St?_l) 



(51) 



where s 



a{3 ,mn 



i(X a -\p) 



Sap.mn, ^a,p being the fields counting the particles in the respective terminals, /(w) 



SmnSa/3 /a(w) is diagonal in both channel (m,n) and terminal (a,/3) indices and describes the energy distribution 
function in the respective terminal. In the simplest situation, when A = and J + = 0, we have four terminals with 
one channel in each of them. The scattering part of the Hamiltonian is Hi = Jj_ (tp^ — ijj)b + J- (ipj — ip s )b, see H28JI (for 
simplicity we ignore the unimportant numerical prefactors). The equations of motion (EoMs) for the participating 
operators read 

id t ips = -id x tp s + J^bS(x) , 
id t ip = —id x ip + J±b5(x) , 

id t b = J x (1>*-il>) + J-(1>l -*!>.)• 

Integrating the first equation over time and then around the point x = we obtain 

i[ip s (0 + )-i> s (0-)) = J_6. 
Acting with idt from the left and using the EoM for the b Majorana one obtains, 

-d t [tlj s (0 + )~t/j s (Q-)] = J-J±(^ - t/j) + J 2 _(^t ~ ^Ps) , 
-d t [yj(0+)~yj(0-)] = J_J ± (^-V s ) + Ji(^ t -V'), 
where the last equation is obtained by symmetry. Now we employ the plain wave decomposition similar to that used 

in EUSEi, 



ctfc for x < 
bk for x > 



ip s (x) 



dk 



Ck for x < 



2ir 1 d k for x > 



Since the dispersion relation of both fcrmion species is trivial, u> = k, we can use u> both for momentum and energy. 
Employing the regularisation scheme ipi = [4>(0 + ) + ^(0 _ )]/2 we obtain 



iuj (du - Cu) = 
-iu) (bu - a u ) = 



,h_J_ 
2 

Ji J_ 



Ji 

2 



(«L. 



Comparing these relations with their adjunct at —to we identify that 



,t 



d) = J — d 



and that 6, 
c. g. 



{duj — Cw) Jj_/J-. Using these expressions we can find both b u and d w as functions of a w and c w , 



1 



" w + i(r x + r_ 

That leads to the following scattering matrix 



(a> + «r_) a w + «r x aL w - iJ±J- c u + iJ±J_ c^_ u 



1 



w + z(r ± + r_) 



uj + iT- iT± —iJ^J- i,J±J- 

iT± u> + ir_ iJ±J- —iJ±J— 

—iJ±J— iJ±J- uj + iT±_ ir_ 

iJxJ- —iJxJ- ir_ a; + ir_i_ 



* a, 1 
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The actual charge transport through the system is conveyed by the charge flavour channel, e. g. by scattering of ip 
fermions across the constriction. The physical picture is similar to that discussed in |2.'"il |: the incoming particles - 
chiral fermions in terminal 1, which are described by ak operators and which have chemical potential \i\ = V - are 
transferred into all other terminals 2-4 (pk, Ck and operators), which are unbiased ^2,3,4 = 0, that is why we have 
to set / = diag(n£, rip, rip, np). Then Ai = A counts particles which leave channel 1. However, the very same fermion 
reappears in the channel 2. Since 1 and 2 are physically one lead we have A2 = —A. We are not interested in change 
of particle numbers in the other channels, that is why A 3 4 = |60l |. Therefore the matrix STs given by 



r_ 



iT x e i2X 
-Uj_J- e lX 



iT x e~ i2X 
to + iT- 

iJj_J-e-' x 
-UiJ-e-' x 



-iJ±J- e- iX iJ ± J^ e- iX \ 

U ± J-e iX -iJj_J-e iX 

lu + iT± ir_ 

ir_ lu + iT± f 



Plugging these relations into (|51fl . folding the integration over energy to the domain [0, 00) and using the properties 
np(-uj) = 1 — np(tu) and n^{— lu) = 1 — n R (cu) immediately leads then to the result (|39[) . 



Appendix C 

Here we establish the relation between our findings and the result of Kindcrmann and Trauzettcl (KT) calculation 
[H. Let us consider Eq. (12) of Ref. [U, 

OO 

InxMX) - ^HVX - T 2 \ 2 ) + T f ^ In {l + T (uj)[f+(1 - D(e lX - 1) + f~(l - f + ){e^ x - 1)]} , 



The /-functions are given by 



j± _ n L/R 



(l-n L/R )e ±lX + n L/R ' 
where u r ^l(lu) = np(uj ± V/2). Obviously 

(l~n L/R )e ±tX 



i-r- 



(1 - n L/R )e ±lX + n L/R 



Taking into account that Tq(lu) = Aj 4 /^ 2 + Aj 2 ), the identification with the KT's impurity strength is 2j = T R - In 
order to proceed we define the object 

ln XifT (A;a) = ^HVA-T 2 A 2 )+T J ^ In {l + aT Q (uj)[f+(l - f~ )(e> x - 1) + /" (1 - f+)(e~ iX - 1)]} , (52) 



where a is a parameter. The derivative of l|52|l with respect to this parameter 

d\n XKT (\;a) _ T °f dcu T (u)[f+(1 - f-)(e iX - 1) + f~ (1 - f+)(e~ iX - 1) 



da J 2tt 1 + oTo(«)[/+(l - /")(e« - 1) + /"(l - f+){e^ x - 1)] ' 



Substituting explicit expressions for the /-functions into the long fraction gives 

T {m)[n L (l - n R ){e- lX - 1) + n R {\ - n L ){e' x - 1)] 

[(1 - n L )e iX + n L ][(l - n L ) e - iX + n R ] - aT Q {u))[n L {l - n R ){e~ iX - 1) + n R (l - n L )(e iX - 1)] 

After simple algebra this simplifies as 

T (u)[n L (l - n R )(e- lX - 1) + n R (l - n L )(e lX - 1)] 
1 + [1 - aT (a;)][ni(l - n R ){e-* x - 1) + n R (l - n L ){e* x - 1)] ' 
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Integrating with respect to a therefore results in 



duj 



\n X K T (\;a) = T I — ln{l + [1 - aT (u;)][n L (l - n R )(e- lX - 1) + n R (l - n L )(e lX - 1)]} + C . 
o 

To fix the constant C, evaluate the above integral at a = 0, see also pfij 

oo 

T / % ln[1 + " n «)( e " A - + ^ - nL ^ glA - !)] = £ HyA - t ' a2) ' 



so that C = 0. The following identity is therefore established (put a = 1): 

oo 

hiXKT(A) = T y|£ln{l + [l- T (w)][ni(l - n R )( e - tX - 1) + » H (1 - n L )(e lA - 1)]} 



(53) 



The two equations 152|) and (|53[l define the same function, which means that the KT statistics, like the g = 1/2 
statistics, is reducible to the generic Levitov-Lesovik formula, its relation to the non-interacting statistics Eq. HOf) 
being 

Xxt(A) = Xo(-A;F;{1-To(w)}) . 
Finally the explicit relation between the KT statistics and our g = 1/2 result is 

X kt(X; V; {ToM» = Xi/ 2 (-A; V/2; {1 - T (u)}) , 
which is a direct consequence of the duality shown in (24|. 



Appendix D 

The calculation starts as usual with the adiabatic potential [sec also Eq. 11321 ] 



d 



U{\±) 



dio r 

2^ 



7+ (sin(A/4)G fc -- + cos(A/4)G^-) + 7 - (cos(A/4)G^ - sin(A/4)G-- 



(54) 



The most compact way to evaluate the inhomogencous GFs entering this expression is through their reduction to GFs 
involving only the resonant level Majoranas a and b. This is accomplished by the following relations: 



G br] 


= H+ 




sin(A_ 


74)% 


-iV sin(A+/4).g+- 


- D bb - 


cos(A_ 


-/4)flf- 


+ D bb + 


cos(A+/4)g+- 




+ i"f- 


~D b - a - 


sin(A_ 


-m- 


-^- a + sin(A + /4) 5 +- 


+ D b ~ 


cos(A_ 


7 4 )%7 


-D ba + 


cos(A+/4).g+- 


G bt; 


= H+ 




sin(A_ 




-D 6 -+ S in(A + /4) 5 +- 


-D bb - 


cos(A_ 


74)%- 


+ D bb + 


cos(A+/4) 5 +- 




+ i"f- 


D ba 


sin(A_ 




-^ 5 - + sin(A + /4) 5 +- 


+ D ba - 


cos(A_ 


74)% 


D ba + 


cos(A + /4).g+- 


c — 


= H+ 




sin(A- 


7 4 )%~ 


- - D ah +sin ( A +/4)fw~ 


-D- b - 


cos(A_ 


74)% 


+ D- b + 


cos(A+/4)%7 




+ ij- 


D-- 


sin(A_ 




- J D-+sin(A + /4) 5 +- 


+ D-- 


cos(A_ 


74)%,- 


- D-+ 


cos(A+/4)%7 




= H+ 


D~ b - 


sin(A_ 


-/4)%- 


-D-+sin(A + /4) 5 +- 


-D- b - 


cos(A_ 


74)%" 


+ D- b + 


cos(A+/4) 5 +- 




+ «7_ 


Daa 


sin(A_ 


74)%- 


- J D-+ S in(A + /4)g+- 


+ D-- 


cos(A_ 


74)%- 


- D-+ 


cos(A + /4).g+- 



Inserting these results into (|54|l leads to 



d 
d>Z 



U{\±) 



1 

~4 
+ 7 
+ 7+7 



dcu c „ 



D 



g vi - cos(A/4)^-+ 5 +- + sin(A/4) J D-+ 5 . 



7)77 



D--g vi - cos(A/4)D-+%- + sin(A/4)D-+ 5 . 



'/'/ 



(D b a D--)g-- sin(A/4)(D b -+ - D"+) fl +- - cos(A/4)(^+ - i%\+)%r] } 



(55) 
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Finally, the calculation of the composite 4x4 matrix object 



V 



D7r D-+ Dr~ D~ + 

bb bb ba ba 

DtT D++ D+- D++ 

bb bb ba ba 

1J ab ^ab ^ aa ^ aa 

D + ~ D ++ D+- D ++ 

^ab ^ab x> aa ^ aa 



can be done by calculation of 



, where = diag(l/o;, — 1/ui, 1/ui, — is the corresponding 



matrix in the absence of the tunnelling couplings and where the corresponding self-energy is given by 



and the components of this object are (we set r± = j±/2 and T± = 7_7 + /2), 



— 



A 2 /u! + iT + (n R + n L - 1) 



-iT + (n L e 



i\/4 



n R e~ 



^^ + {l-n R )e a ' A + {l-n L )e- fx ' i -A 2 /uj + iT + (n R + n L -l) 



n R — til 7iLe lX / 4 - 

(1 - n R y x / A - (1 - n L )e~* X/4 n R 



n R e 

- n L 



-iA/4 



n L - n R 
-(1 - n R )e tX / 4 + (1 - n L )e- lX / 4 



-n L e fx / A + n R e~ lX / A 
n L - n R 



A 2 /u + iT-(n R + n L - 1) 



-ir_(n L e a / 4 + n fl e-' a / 4 ) 



ir_ (1 - n R )e lX / 4 + (1 - n L )er lX l A -A 2 /uj + iT-(n R + n L - 1) 



Unsurprisingly, an inversion of the matrix (2? < -°)) -1 — E g and substitution of the result into l|55|l leads to exactly the 
same expression for the statistics l|49|) up to the transmission coefficient, which is now more involved and which is the 

same as calculated previously in the context of the non-linear I — V of the same set-up I23 ■ 
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